A Correlation between the Emission Intensity of Self- Assembled Germanium Islands 
and the Quality Factor of Silicon Photonic Crystal Nanocavities 



N. Hauke/'QS. Lichtmannecker, 1 T. Zabel, 1 F. P. Laussy, 1 A. Laucht, 1 
M. Kaniber, 1 D. Bougeard, 2 G. Abstreiter, 1 ' 3 J. J. Finley, 1 and Y. Arakawa 3 ' 4 

1 Walter Schottky Institut and Physik Department, 
Technische Universitat Miinchen, Am Coulomb-wall 4, D-85748 Garching, Germany 
2 Institut fur Experimentelle und Angewandte Physik, 
Universitat Regensburg, D-93040 Regensburg, Germany 
^Institute for Advanced Study, TU Miinchen, D-85748 Garching, Germany 
4 Institute of Industrial Science, Institute for Nano Quantum Information Electronics, 
The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 153-8505, Japan 
(Dated: January 15, 2013) 

We present a comparative micro-photoluminescence study of the emission intensity of self- 
assembled germanium islands coupled to the resonator mode of two-dimensional silicon photonic 
crystal defect nanocavities. The emission intensity is investigated for cavity modes of L3 and 
Hexapole cavities with different cavity quality factors. For each of these cavities many nominally 
identical samples are probed to obtain reliable statistics. As the quality factor increases we observe 
a clear decrease in the average mode emission intensity recorded under comparable optical pumping 
conditions. This clear experimentally observed trend is compared with simulations based on a dissi- 
pative master equation approach that describes a cavity weakly coupled to an ensemble of emitters. 
We obtain evidence that reabsorption of photons emitted into the cavity mode is responsible for 
the observed trend. In combination with the observation of cavity linewidth broadening in power 
dependent measurements, we conclude that free carrier absorption is the limiting effect for the cavity 
mediated light enhancement under conditions of strong pumping. 

PACS numbers: 42.60.Da 42.70.Qs 78.55.-m 42.50.Ct 



I. INTRODUCTION 

The development of an efficient silicon (Si) based 
light source is of great interest in the information tech- 
nology industry since it would pave the way towards 
optical interconnects with CMOS compatibility. If 
achieved, this major goal would dramatically enhance 
the signal processing speeds currently achievable in Si- 
microelectronicsP^ Due to its indirect bandgap light 
emission in crystalline Si requires the participation of 
phonons to conserve crystal momentum. Typically, this 
leads to a very low internal quantum efficiency and Si is, 
therefore, rarely used as active light emitting material. 6 
One approach that has been explored to realize an ef- 
ficient Si based light source is to enhance the material 
radiative emission efficiency by exploiting cavity quan- 
tum electrodynamic effects using photonic crystal (PhC) 
nanocavities. Enhanced photoluminescence (PL) has 
been recently reported in crystalline Si PhCs-^ as we ll 
as for germanium islands (Ge-islands) embedded in a Si 
PhC nanocavityP^HH The future development of efficient 
Si-based light sources and, potentially, even a CMOS 
compatible laser would revolutionize information tech- 
nologies. However, detailed investigations of the nature 
of the light matter coupling have not been performed, 
either theoretically or experimentally. Therefore, little 
is known about the emissive properties of Ge-islands in 
nanocavities. Furthermore, PhC nanocavities might pro- 
vide a route to enhance optical activity to a level where 
optical properties of single Ge-islands can be investi- 



gated. 

In this article we report on the investigation of the 
cavity enhanced emission from self-assembled Ge-islands, 
which are grown by molecular beam epitaxy (MBE) and 
are embedded in PhC nanocavities. We begin by compar- 
ing the PL intensities from L3 and Hexapole PhC cavity 
modes under comparable conditions of optical pumping 
as a function of their quality (Q)-factors. Averaged over a 
large number of different and nominally comparable cavi- 
ties, a very general finding of our work is that an increase 
of the average mode PL-intensity is observed as the mode 
Q-factor becomes smaller. The average emission inten- 
sity saturates for Q-factors less than s» 600. In order 
to understand this very clear experimental observation 
we introduce a cavity quantum electrodynamics (QED) 
model based on a dissipative master equation for an en- 
semble of emitters, where each emitter is in a highly ex- 
cited state to account for the strong optical pumping and 
weakly coupled to the cavity mode to account for the low 
oscillator strength. The model exhibits various regimes 
with strongly different dependencies between mode in- 
tensity and Q-factor, depending on the spectral emit- 
ter ensemble distribution and the emitter-cavity detun- 
ing. By comparing this model with the experimentally 
observed trend we obtain evidence that reabsorption of 
photons emitted into the cavity mode is responsible for 
the observed Q-factor dependence of the intensity - an 
effect that ultimately limits the radiative efficiency. By 
fitting the predictions of our model to our data we ex- 
tract a photon reabsorption time of r a b s = 2.3 ± 0.4 ps 



in very good agreement with the photon reabsorption 
time r* bs = 1.9 ± 0.3 ps extracted from power dependent 
measurements of the cavity mode linewidth. This sug- 
gests that free carrier absorption (FCA) takes place un- 
der the conditions of strong optical pumping in the cavi- 
ties. Hence, cavity modes with high Q-factors and, thus, 
long photon lifetimes exhibit a reduced internal radiative 
efficiency under conditions of strong pumping compared 
to cavity modes with low Q-factors. 



II. SAMPLE STRUCTURE AND EXPERIMENT 

The sample consists of a 2D PhC that is fabricated into 
a freely suspended Si slab, which contains a single layer 
of self-assembled MBE grown Ge-islands as illustrated 
in fig. [l](a). The fabrication starts from a silicon-on- 
insulator (SOI) wafer provided by Soitec^with a 220 nm 
thick crystalline Si layer on top of a 3 pm thick layer of 
buried SiC>2. Before growth, the crystalline silicon layer 
is thinned to 50 nm using isotropic wet chemical etching 
with a mixture of 60% nitric acid (HNO3) and 0.04% hy- 
drofluoric (HF) acid. After the transfer into the MBE the 
native oxide is thermally removed by heating to 760° C 
for several minutes. The crystal growth is then initialized 
with a 85 nm Si buffer layer grown at 520-560°C. Follow- 
ing this, six monolayers (ML) of Ge and a 135 nm thick 
Si capping layer are deposited at temperatures of 430° C 
and 410°C, respectively. For structural investigations an 
uncapped layer of islands was grown on the sample sur- 
face with the same growth conditions used for the capped 
island layer. Fig. [I] (b) shows an atomic force microscope 
(AFM) image of the surface islands, revealing a bimodal 
island distribution with smaller "pyramids"^ and larger 
"domes'^. The emission from the two types of nanos- 
tructures is spectrally distinct!^ and, for the PhC cavities 
fabricated in this work, the cavity modes are tuned into 
resonance with the dome emission at 0.92 eV. Thus we 
focus on the interaction between the Ge-dome-islands, 
which we refer to in the following as "Ge-islands", and 
the PhC nanocavity modes. 

After growth, photonic crystal nanostructures were re- 
alized using electron beam lithography and subsequent 
SF6/C4F8 reactive ion etching (RIE) to define hexago- 
nal lattices of air holes with three different periods of 
a± = 330 nm, 02 = 360 nm and 03 = 390 nm. The scan- 
ning electron microscope (SEM) image in fig. [I] (c) shows 
a typical PhC with a lattice constant of 03 = 390 nm con- 
taining a L3 PhC nanocavity.^ As a final processing step 
the underlying Si02 is selectively removed by HF acid to 
form a freestanding slab membrane. 

Optical measurements were performed using a micro- 
photoluminescence (/iPL) spectroscopy setup. The sam- 
ple was placed in a liquid helium flow-cryostat for low 
temperature investigations. To excite the sample we used 
a continuous-wave (cw) diode-pumped solid-state-laser 
emitting at XLaser = 532 nm, which is focussed by a lOOx 
microscope objective (NA=0.5) to a spot size with a di- 
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FIG. 1. (Color online) (a) Left panel: Schematic cross- 
sectional representation of the photonic crystal nanocavity 
structures investigated. Right panel: Layer sequence in the 
active region, (b) Atomic force microscope image of the Ge- 
islands investigated, (c) SEM image showing a L3 photonic 
crystal cavity from the top. (d) /iPL spectrum from Ge- 
islands in the unpatterned region at T = 25 K. The red line 
is a Gaussian fit to the data. 



ameter of rs 0.8 pm. The resulting PL signal is collected 
through the same objective and dispersed by a 0.32 m 
imaging monochromator equipped with a 600 lines/mm 
grating and a liquid nitrogen-cooled InGaAs linear diode- 
array. 

In Fig. [I] (d) we present a typical //PL spectrum of the 
Ge-islands emitting in the unpatterned region of the sam- 
ple. The data, obtained at T — 25 K, can be fitted well 
by a Gaussian peak with a full width at half maximum 
(FWHM) of s» 50 meV as indicated by the red line. Due 
to carrier diffusion our excitation spot size of sa 0.8 pm 
leads to a region with a FWHM of ~ 1.5 pm, that gener- 
ates PL-signal, as obtained by performing /iPL with spa- 
tially separated excitation and detection spots. By com- 
paring this finding with AFM measurements performed 
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on uncapped surface islands, as shown in fig. [T] (b), we es- 
timate that w 100 Ge-islands are optically excited by our 
laser. Hence, the observed spectrum represents the sum 
of the individual spectra of N « 100 single Ge-islands. 

In the following, we will compare the emission inten- 
sity of cavity modes with various Q-factors to study the 
cavity-emitter coupling. For all experiments reported in 
this paper the nominal emission intensity was recorded 
using an optical excitation power density of 600 kW/cm 2 
in order to allow a comparative study. Furthermore, it 
is required that the modes couple spectrally to the same 
type of islands and, thus, are in a specific spectral win- 
dow. To do this we tuned the lattice constant of our 
PhC structures during the fabrication process to coarsly 
adjust the cavity modes and varied the air hole diameter 
to fine tune the energy. In fig. [5] (a) we plot ^PL spectra 
recorded at T = 25 K from a series of L3 PhC nanocav- 
ities with a lattice period of ci2 = 360 nm and different 
air-hole radii, increasing from bottom to top. Six distinct 
emission lines from the L3 cavity can be observed, five 
of which are clearly observable in the figure (Ml is the 
fundamental cavity mode, M2-M5 are the higher energy 
modes) . When increasing the air hole diameter the mode 
emission shifts systematically to higher energies. For our 
analysis we will only consider mode emission recorded 
from samples in the spectral region between 0.915 eV 
and 0.935 eV, highlighted by the yellow shaded region 
in fig. [2] (a). As can be seen in the figure, we tune the 
emission of M2 to M5 into this spectral reference region 
via the fine tuning method alluded to above. For spec- 
tral coarse tuning of the mode emission we use the PhCs 
with lattice constants of ai = 330 nm and 03 = 390 nm 
to bring Ml and M6 into the spectral region of interest, 
respectively. 

In order to quantitatively compare the intensity of 
the PL emission we have to take the mode volume of 
the nanocavity modes into account. This quantity in- 
fluences the coupling strength between emitter and cav- 
ity and, thus, the cavity mode /aPL intensity. In or- 
der to exclude this effect we consider only cavity modes 
with similar values of the mode volume in our anal- 
ysis. Therefore, we performed numerical finite differ- 
ence time domain (FDTD) simulations^ which show that 
Ml, M2, M3 and M6 have comparable mode volumes 
of Vmodo = 0.65 ± 0.09 (A/n) 3 , whilst M4 and M5 ex- 
hibit significantly lower mode volumes of 0.30 (X/n) 3 and 
0.39 (X/n) 3 , respectively. Therefore, we excluded M4 and 
M5 from our analysis and consider only Ml, M2, M3 and 
M6. In addition to the L3 cavities, we also fabricated 
hexapole cavities 20 emitting in the spectral reference re- 
gion and showing a mode volume of V mo dc = 0.63 (X/n) 3 
for the fundamental dipole mode (Ml hex) of this struc- 
ture. Hence, we include this mode in the evaluation pre- 
sented below. In fig. [2] (b) we plot representative mode 
emission spectra in the spectral region of interest from 
0.915 eV to 0.935 eV for a number of different cavity 
modes: (i) the fundamental L3 cavity mode (Ml), (ii) 
the first higher energy L3 mode (M2), (hi) the second 
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FIG. 2. (Color online) (a) fj,PL spectra recorded at T = 25 K 
from a series of L3 PhC cavities with different air-hole radii, 
increasing from bottom to top. (b) Mode emission in the spec- 
tral range from 0.915 eV to 0.935 eV [marked yellow in (a)] for 
different cavity modes: (i) fundamental L3 cavity mode (Ml), 
(ii) first higher energy L3 mode (M2), (iii) second higher en- 
ergy L3 mode (M3), (iv) fifth higher energy L3 mode (M6) 
and (v) fundamental dipole mode of a hexapole cavity (Hex 
Ml). 



higher energy L3 mode (M3) , (iv) the fifth higher energy 
L3 mode (M6) and (v) the fundamental dipole mode of a 
hexapole cavity (Hex Ml). The PL intensity scale is the 
same for all graphs. 

In addition, we need to take the far field emission pro- 
file of the different cavity modes into account, since the 
far field profile influences the collection efficiency ?7 co ii of 
our optical detection system. As a result, the fraction 
of light emitted to one hemisphere that is collected by 
our microscope objective differs for different cavity mode 
profiles. We obtained the values of n co \\ for the differ- 
ent cavity modes investigated using FDTD simulations. 
Table U summarizes the simulation results for the cavity 
modes we include in our analysis. Here, a is the lat- 
tice constant used for coarse shifting the mode emission, 
Kuode the mode volume and r) co u the photon collection 
efficiency. 

Finally, we need to exclude measurement errors such 
as imprecise positioning of the excitation spot on the 
cavity or varying fabrication quality of different PhCs, 
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TABLE I. Overview of the simulated properties of the PhC 
cavity modes included in our mode intensity versus Q analy- 
sis. 



Mode 


a(nm) 


V mo dc(A/n) 3 


??coll 


Ml 


330 


0.74 


0.197 


M2 


360 


0.73 


0.108 


M3 


360 


0.57 


0.437 


M6 


390 


0.63 


0.468 


Ml hex 


330 


0.63 


0.466 



(Q < 600) are similar. In contrast, for j a < 1.5 meV 
(Q > 600), we observe a progressive and systematic de- 
crease in J mo dc with decreasing cavity linewidth (increas- 
ing Q-factors). For the fundamental mode of the L3 cav- 
ity (Ml), which is the mode with the smallest linewidth, 
^modc is strongly reduced to approximately 1/5 of the 
intensity of Hex Ml and M3. In order to understand 
this very clear experimental observation, we introduce in 
the following section a cavity quantum electrodynamics 
(QED) model based on a dissipative master equation. 
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FIG. 3. (Color online) Average mode emission intensity as a 
function of inverse cavity Q-factor (cavity linewidth 7 a ). Each 
point is obtained by averaging over several modes emitting in 
the spectral range between 0.915 eV and 0.935 eV. The error- 
bars display the standard-deviation. 



all of which would influence the intensity of the PL sig- 
nal. Hence, we measured 3-7 different photonic crystal 
nanocavities for each of the mode-types emitting in the 
spectral window from 0.915 eV to 0.935 eV. The mode 
intensity is extracted by fitting a Lorentzian peak to the 
spectral profile and averaging the peak areas for modes 
of the same type. The background stemming from the 
uncoupled Ge-island emission or from spectrally closely 
spaced modes (M3 and M4) is subtracted to separate the 
cavity emission from the background. 

In order to obtain the emitted mode intensity we cor- 
rect the measured intensity for the mode volume (Vmode) 
related change of the number of islands coupled to the 
cavity, for the angular collection efficiency (recoil) and 
for the mode degeneracy Dp3 In fig. [3] (a) we plot the 
obtained average emitted mode intensities (I mo dc) as 
a function of inverse Q and cavity mode linewidth 7 a 
(Q = uJcav/la)- The error bars account for the stan- 
dard deviation of the intensity and Q-factor distribution, 
since each point is the average of several independent 
measurements performed on a number of cavities. The 
average emission intensity of cavities with 7 a > 1.5 meV 



III. THEORY 

In this section we develop a cavity-QED model that de- 
scribes our system of an ensemble of Ge-islands coupled 
to a cavity mode. Ge-islands are expected to have a very 
weak transition dipole moment, due to the spatial sepa- 
ration between electrons captured at the Si-Ge interface 
and the holes localized in the Ge-islands (type II band 
alignment).^ Hence, the coupling parameter g that de- 
scribes the interaction between emitter and cavity, is ex- 
pected to be very small. In the following the parameters 
will be denoted as illustrated schematically in fig. [4] (a) , 
where N is the number of Ge-islands coupled to the cav- 
ity mode, 76 is the FWHM of the emission of a single 
Ge-island, a is the FWHM of the spectral island ensem- 
ble distribution, 7 a is the FWHM of the cavity emission 
and A is the effective spectral detuning between the cen- 
ter of the cavity and the center of the island ensemble. To 
capture the essential physics that determines the dynam- 
ics of such a system, i.e. N Ge-islands weakly coupled to 
a cavity mode, we turn to the simplest possible picture 
provided by a dissipative master equation ifid t p = Cp 
for the density matrix p with Cp = [H, p] + C 1 p + Cpp, 
where H is the Hamiltonian that describes the dynamics 
of the islands-nanocavity system: 



N 



N 



H = H h . cc + J2^9t(b]a + b^ 



(1) 



Hf ree = hui a a>a + Y^iLi huib\bi being the free dynamics 
of the modes and 

£ 7c p = -~ (2cpc f -Sep- pc>c) , (2) 

is the Lindblad operator for the cavity modes, where 
c = a, bi is the decay of the modes. Pumping is included 
in the same way, as an incoherent source of excitation 
with Lindblad term Cp^p = ^f(2b\pbi — bib\p — pbib\). 
The nonlinear term b\b\bibi describes phase space filling, 
by approximating the real energy level structures of the 
island^ to be an an equally spaced ladder of levels^ with 
spacing Vt for the i th dot. A mean field approximation 
{b\bia*bi) s» (6j6i)(a 1 '6i) is performed that allows us to 
truncate the equations of motion self-consistently. The 
main effect of this term is to provide an effective detun- 
ing between the cavity mode and the spectral center of 
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the emitter ensemble distribution, as well as an effective 
broadening of the emitter line. The steady state popu- 
lation in the cavity n a = lim^oo Tr(p(i)a^a) is then ob- 
tained, from which follows the number of photons emitted 
per unit time J mo de = laJia- This is the main quantity 
of interest in our experiment. A closed form expression 
can be obtained for a distribution of islands (a Gaussian 
distribution is non-integrable) that shows the combined 
effect of detuning A = uj a — X^iLi u i between the cavity 
mode and the average position of the islands and their 
distribution cr: 

j _ p b la 4 9c 2 ff 

-'mode ^ / . -pi , _\ ■ \°) 

lb (7a + T 6 + a) 7 a 

Here, we have retained the leading term only in the cou- 
pling strength <?, since it is very small in our systems, and 
have defined the effective coupling: 

a Ng 2 

9eS = ~. T2 > ( 4 ) 



where A = \J A 2 + [(h^ b%Vi)/N} 2 is the effective de- 
tuning that includes the interactions, and is the effec- 
tive broadening of a single emitter that, following Bose 
statistics, reads 7b — Pb- These two new parameters are 
however to be considered the natural and independent 
ones that describe the system, rather than the micro- 
scopic ones from which they stem. A and 7b do not play 
a direct role anymore. 

In terms of these effective parameters, eqn. Q is 
closely related to that of a single emitter coupled to a cav- 
ity (and reduces to it when N = 1, Vi = and a = Ojp^. 
The first term in eqn. Pb/Tb, is the effective pop- 
ulation that builds up among the islands. The second 
term, ^ a /{la + + cr), is the fraction of these excita- 
tions which is available to excite the cavity. The third 
term, 4g 2 ff /7 a , dominates the dynamics and governs the 
intensity-linewidth trend. 

In fig. [4] (b) we plot eqn.Q for various limiting cases of 
the cavity QED light-matter coupling: (i) Cavity at reso- 
nance (A = 0) and (ii) at finite detuning (A = 1.5 meV) 
to a spectrally narrow emitter ensemble with an effec- 
tive ensemble linewidth of Tb + cr = 0.01 meV. In 
curve (iii) we plot eqn. ^ for a cavity at large detun- 
ing (A = 30 meV) to a broad emitter ensemble distri- 
bution with Tb + cr = 50 meV. These parameters corre- 
spond closely to the expected reality for our system cf. 
fig. [11(d)]. For very weak coupling (g 0), the inten- 
sity of the emitted light increases with the quality of the 
cavity ("f a — > 0) if the cavity mode is placed at resonance 
to a spectrally narrow emitter ensemble (A = 0). This 
can clearly be seen by curve (i) in fig. [4] (b) . At non-zero 
detuning with finite spectral mismatch A the effective 
coupling becomes larger with decreasing Q (increasing 
7 a ) due to an increasing overlap between the detuned 
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FIG. 4. (Color online) (a) Schematic illustration of the pa- 
rameters used for the theoretical model of N emitters coupled 
to a cavity: Cavity linewidth j a , emitter linewidth 75, emitter 
distribution a and effective ensemble cavity detuning A. (b) 
Theoretical trends in various limiting cases of the cavity QED 
light-matter coupling: (i) Cavity at resonance (A = 0) and (ii) 
at finite detuning (A — 1.5 meV) to a spectrally narrow emit- 
ter ensemble with a FWHM of T b +a = 0.01 meV. (iii) Cavity 
at large detuning (A = 30 meV) to a broad emitter ensem- 
ble with an effective ensemble linewidth of Tt + a = 50 meV, 
which corresponds to our system. 



cavity mode and the collective set of islands, with effec- 
tive linewidth Tb + cr. Hence, making the cavity Q worse 
gives rise to an increase of the mode emission intensity. 
For very low Qs (large 7 a ) the intensity decreases since 
the coupling between the emitter ensemble and the cav- 
ity gets weak. This is shown representatively by curve 
(ii) in fig. [4] (b). The trends shown in curves (i) and 
(ii) can be observed only when the effective ensemble 
linewidth Tb + cr of the island-ensemble is comparable 
to the linewidth 7 a of the cavity mode. In strong con- 
trast, our system has an ensemble linewidth of « 50 meV, 
which is more than an order of magnitude larger than the 
cavity linewidth. For this case our model predicts that 
the emitted intensity should exhibit a plateau for the 
regime of cavity mode linewidths we observe in our ex- 
periment (0.3 meV < j a < 3 meV), quite independently 
of detuning. This is shown by curve (iii) in fig. [4] (b). 
Clearly, experimentally we do not observe the theoreti- 
cally predicted plateau for the high-Q region, but rather 
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a decrease in emission intensity as presented in fig. [3p3 
So far, our cavity-QED model takes only photon emission 
from the islands into the cavity mode into account, while 
in the real experiment we have additional effects, such as 
photon reabsorption, annihilating the photon before it es- 
capes the cavity. With increasing cavity Q the time the 
photon remains inside the cavity increases and, hence, 
absorption effects are expected to play an increasingly 
important role for high Q cavities. This is supported by 
the fact that Ml is expected to have a high Q-factor of 
Qsim ~ 70 000 according to our FDTD simulations, but 
we observe a significantly smaller Q-factor in our exper- 
iment of Q oxp ps 2 800. 

Due to the spatial separation between electron- and 
hole-wavefunction and the indirect optical transition in 
/c-space we expect the resonant absorption of the Ge- 
islands to be too weak to explain the Q-factor saturation 
in our experiment. This expectation is supported by the 
observation that the Q-factor of the nanocavity emission 
is independent of the spectral position relative to the Ge- 
island ensemble. If reabsorption of photons by Ge-islands 
would play a dominant role, we would expect to observe 
an increase of Q-factors for cavity modes emitting at the 
low energy side of the ensemble, since the probability 
of photon reabsorption decreases there. However, free 
carrier absorption (FCA) can cause photon reabsorption 
since we generate a large density of charge carriers in 
the vicinity of the nanocavity, as we optically excite with 
a high power density of 600 kW/cm 2 . Such a high op- 
tical pumping intensity is required to observe a strong 
PL signal. We take reabsorption of photons into account 
as a competing process between the escape of a photon 
out of a cavity with a photon escape time r esc , given by 
the intrinsic Q-factor of the cavity, and the photon reab- 
sorption time T a bs- The total photon loss time r a , given 
by l/r a = 1/t csc + l/r a b s , determines the experimentally 
observed Q-factor and, therefore, the cavity linewidth j a 
measured in PL. Hence, the emitted mode intensity i^Jode 
in the presence of photon reabsorption can be expressed 



^modc(7a) — -^modc (7csc ) ' 



l/r« 



l/ T csc + 1/ T abs 
Tabs 



^mode(Ta Tabs) (1 )■ 
la 



(5) 



-'mode (7<») denotes the emitted PL intensity of a cav- 
ity with a linewidth of j a when including absorption, 
^mode(7esc) denotes the emitted PL intensity of the same 
cavity without absorption and, hence, a smaller cavity 
linewidth 7 CSC which is solely determined by l/r osc . The 
last term describes the competing process between the 
escape of a photon out of the cavity and photon reab- 
sorption, where 7 a b s is the cavity linewidth given by the 
absorption. 



IV. DISCUSSION 

In fig. [5] (a) we present the theoretically predicted 
trend of the mode emission intensity as a function of cav- 
ity linewidth j a for our system in the absence of photon 
reabsorption [dotted line, cf. fig. [4] (b)]. Furthermore, 
we plot a fit (full line) of the model to the experimental 
datapoints (grey, points) from fig. [3j using the same pa- 
rameters and now including reabsorption of photons as 
described by eqn. [5j There, we used 7 a b s and the normal- 
ization of I mo dc as fitting parameters. The theoretically 
predicted trend follows our experimental data by showing 
a decrease of emission intensity with decreasing j a (in- 
creasing Q-factors). Obviously, we do not expect to ob- 
serve mode emission for linewidths "f a smaller than 7 a b SJ 
as highlighted by the red area. From the fit we extract 
7abs = 0.27 ± 0.05 meV, which corresponds to a photon 
reabsorption time of r a b s = 2.3 ± 0.4 ps. The presence of 
FCA is supported by the fact that we observe, with in- 
creasing optical pumping power, a broadening of the cav- 
ity modes. This is an effect, which El Kurdi et a/P^have 
observed at similar excitation power densities for Si pho- 
tonic crystal nanocavities with embedded Ge-islands. An 
example of such pump power dependence investigations 
for M2 is shown in fig. [5](b), where we plot the mode 
linewidth of M2, 7 ^ 2 , as a function of optical excitation 
power density. In our system, at low optical excitation 
power densities, the linewidth of M2 stays constant at 



7, 



M2 



0.74 ± 0.02 meV (marked by the horizontal dot- 
ted line) , as it is solely determined by the intrinsic photon 
escape time r esc of the cavity mode. When increasing the 
excitation power density above 250 kW/cm 2 , we observe 
a strong increase in linewidth. At 600 kW/cm 2 (marked 
by the green vertical line), which corresponds to the op- 
tical excitation power density we used in our compar- 
ative experiment, 7^ 2 is increased to 0.97 ± 0.02 meV. 
When further increasing the optical excitation power, 
the absorption coefficient is increased due to heating 
caused by the excitation laser. This leads to a strong 
increase of optically generated charge carriers inside the 
cavity and, as a result, 7™ increases strongly. Since 
we can extract the intrinsic linewidth of M2 from the 
measurements at low power, we can extract the reab- 
sorption 7 abs , using 7 abs = 7^ 2 - 7^ 2 . In fig. [5] (c) we 
plot 7 a b s as a function of optical excitation power den- 
sity. At 600 kW/cm 2 (marked by the green vertical line) , 
which corresponds to the optical excitation power den- 
sity we used in our comparative experiment, we obtain 
7abs = 0-23 ± 0.04 meV, corresponding to a photon reab- 
sorption time of r* bs = 1.9 ± 0.3 ps. This is in remarkable 
agreement with r a b s = 2.3 ± 0.4 ps, as extracted from 
the fit in fig. [5] (a). Thus, these observations lend strong 
support to the argument that FCA limits the emission 
enhancement. A major conclusion of our work is that 
the overall emitted intensity from optically pumped PhC 
nanocavities with embedded Ge-islands is mainly limited 
by FCA, caused by the high charge carrier density from 
strong optical excitation. Concerning the external quan- 
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FIG. 5. (Color online) (a) Mode emission intensity as a func- 
tion of cavity linewidth y a '- (dotted line) Theoretically pre- 
dicted trend for our system with an emitter ensemble distribu- 
tion with Tb + cr = 50 meV and an ensemble cavity detuning of 
A = 30 meV, in the absence of photon reabsorption and (solid 
line) including reabsorption of photons with an absorption of 
Tabs = 0.27 ± 0.05 meV together with the experimental dat- 
apoints from fig. [3] (b) Cavity mode linewidth of M2, y^ 12 , 
as a function of optical excitation power density. The vertical 
line marks the optical excitation power density we used in our 
comparative study (600 kW/cm 2 ). The horizontal dotted line 
marks the intrinsic linewidth 7^,? of M2. (c) 7 a b s , extracted 
from 7^ 2 , as a function of optical excitation power density. 
The vertical line marks the optical excitation power density 
we used in our comparative study (600 kW/cm 2 ). The hor- 
izontal dotted line marks the extracted absorption linewidth 
°f 7abs = 0.23 ± 0.04 meV at this power. 



turn efficiency, high-Q cavities are actually performing 
worse than low-Q cavities due to the increased probabil- 
ity of reabsorption of photons emitted into the cavity. 



CONCLUSION 



We reported a comparative study of PL emission from 
PhC nanocavities with embedded Ge-islands by low tem- 
perature fiPh spectroscopy. First we investigated a num- 
ber of different L3 and Hexapole PhC cavity modes with 
various Q-factors. For a valid comparison we considered 
only mode emission in a well defined spectral region, 
corrected for the mode volume related change of emit- 
ters pumping the cavity mode and the far field radiation 
pattern. With increasing Q-factors we observed a de- 
crease of the PL-intensity emitted by the cavity mode. 
We then introduced a cavity-QED model based on a 
dissipative master equation to understand the dynamics 
of an ensemble of emitters, which are in highly excited 
states due to strong optical pumping and in very weak- 
coupling with the cavity mode. With this model we could 
identify various regimes of PL intensity versus Q-factor 
trends, depending on the emitter ensemble distribution 
and the spectral emitter cavity detuning. By compar- 
ing the theoretically predicted trend of our system with 
the experimental data, we concluded that reabsorption of 
photons emitted into the cavity limits the emission en- 
hancement via the cavity mode. We extracted a photon 
reabsorption time of r abs = 2.3 ± 0.4 ps and concluded, 
by comparing this value to the photon reabsorption time 
r abs = i 0-3 ps extracted from power dependent mea- 
surements, that the absorption is caused by free carri- 
ers. Therefore, cavity modes with high Q-factors and, 
hence, long photon lifetimes exhibit a reduced radiative 
quantum efficiency compared to cavity modes with low 
Q-factors. 
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